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N omenclature 


Natural units used in this work are such that H = c = 1. 
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Lorentz factor, , 
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Abstract 


String fragmentation models such as the Lund Model fit experimental particle 
production cross sections very well in the high-energy limit. This paper gives 
an introduction of the massless relativistic string in the Lund Model and shows 
how it can be modified with, a simple assumption to produce formulas for meson 
production cross sections for space radiation research. The results of the string 
model are compared with inclusive pion production data from, proton-proton 
collision experiments. 

1. Introduction 

The Lund model is a (1 + 1) massless relativistic string fragmentation model which is 
modified with a simple assumption in this paper to produce formulas for meson particle 
production cross sections for HZETRN. The idea of modelling hadronic systems with 
strings went back to the 1960s (ref. 1). The general skepticism on the extra dimensions 
predicted by string theory and the advent of QCD in the 1970s put string theory out of 
commission for the next decade. Interests in string theory rekindled in the 1980s when 
Green and Shwarz showed that string theory is anomaly free and probably finite to all 
orders in perturbation theory. Today many advances in string theory have been made 
on the theoretical front. Despite all the excitement string theory has generated in the 
high-energy theory community, no evidence exists that string theory will make any low- 
energy limit prediction to be tested by experiments in the near future. Supersymmetry 
and superstring may be testable by experiments in the TeV scale but perhaps not directly 
relevant to the focus of space radiation research in the GeV scale. The interest lies 
mostly in the low-energy limit of a nonsupersymmetric string theory, with the possible 
exception of understanding the Greisen-Zatsepin-Kuzmin (GZK) cutoff in the ultrahigh- 
energy cosmic ray spectrum using TeV strings (ref. 2). 

String fragmentation has been a work horse in analyzing high-energy particle pro- 
duction. Isgur proposed a fluxtube model in which the color force field is thought as 
a string-like object (refs. 3 and 4). B. Andersson implemented the idea of string frag- 
mentation formally in the Lund Model (refs. 5, 6, and 7), which is implemented in the 
JETSET and PYTHIA Monte-Carlo programs (ref. 8). If a Monte-Carlo cross-section 
program is put into a target code, it will not run in a short time. Therefore we need 
simple parameterizations of cross-section formulas, which is the aim of the present work. 
In this paper, we first review the original Lund Model and expand some of the deriva- 
tions of reference 5. Later we show how to modify the Lund Model by inserting a simple 
assumption to generate formulas for meson production cross sections. 

Appendix A defines the basic kinematic notations. Appendix B explains the basic 
concepts of the Lund Model and derives the invariant amplitude formulas. The formulas 



ill the appendixes are taken from reference 5, chapters 6 to 8. This work merely serves 
to focus on a subsection of reference 5 that is relevant to our discussion and expands the 
derivations. Section 2 of this paper explains a new idea on how to use a simple ansatz to 
obtain formulas for production cross sections in the Lund Model. 

2. Cross-Section Formulas 


The basic result of the Lund Model is the “Area Law” which is summarized as (ref. 5) 
dPext = ds ^ (1 - z) a e~ br 

Poj\ e~ bArest 

and derived in appendix B. The symbols T and A rest are Lorentz invariant kinematic 
variables. The variable T defines the surface of constant proper time along which the string 
is broken. Traditionally the Lund Model results in equations (1) and (2) are incorporated 
into Monte Carlo simulation programs such as JETSET and PYTHIA to compute cross 
sections. Since the goal of the HZETRN program is to derive simple formulas for the cross 
sections, numerical calculations are not wanted. In this section, some simple assumptions 
are made to derive analytical results. 
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Equation (2) resembles the quantum mechanical result 
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In principle, A tota i = T + A rest can also be used in equation (5) in place of A rest . The only 
difference is a proportionality constant exp(— bT) in the invariant amplitude |.M| 2 . Hence 
the assumption in equation (5) is the simplest. The string breaks up along a surface of 
constant proper time T which is determined by the distribution (ref. 5) 

H(T) =CT a e~ br . (6) 


By setting dH/dT = 0, we can easily show that the maximum of H(T) is (ref. 5) 


r 


max 


a 

b 


(7) 


2 



The mean proper time is (ref. 5) 


J 0 °°r//(r)dr a + i 

Jo” H(T)dr ~ b 


( 8 ) 


The fact that a and b are constants in T as shown in appendix B does not exclude the 
possibility that they are functions of rapidity. Since string fragmentation is a stochastic 
process that occurs over a surface of constant proper time on the average, < T > is 
expected to be an absolute constant; this means that 


a = a 0 f(y) - 1 (9) 

b = 60/(3/ ) (10) 


for some function of rapidity f(y). The goal of a consistent string theory is to calculate 
a 0l b 0 , and f(y) analytically. In this paper, the simple assumption is that 


f(v) 


2 niT sinh y 2 p z 



x, 


( 11 ) 


where nir = \Jmft + p^. The kinematic identities of equation (11) are taken from ref- 
erence 9. The ansatz f(y) = x is motivated by a posteriori considerations when fitting 
experimental data. In NASA space radiation research, cosmic particles are highly en- 
ergetic with production mostly forward. In the center-of-mass frame, the longitudinal 
momentum is p z and the total energy of the colliding particles is y / s'. According to the 
Fermi Golden Rule, 

Together with equations (10) and (11), the differential cross section is found to be 

3— = Ae~ Bx (13) 

ax 

for some constants A and B. This prediction is consistent with experimental data from 
references 10 and 11 as shown in figures 1 and 2. 


The spectral distribution can be obtained from the cross-section formula in reference 12 

Her r^max eft CT 

lb '"-''l deE W 3in6 (14) 

By using approximations p z « E and x p z /\fs in the high-energy limit and equa- 
tion (12) in conjunction with equation (14), the spectral distribution is found to be 

2tt(1 - cos 0 max ) p z \M\ 2 
cEe- kE 



3 


(15) 



where c and k are constants. Unfortunately no experimental data exist for do/dE for 
pion production in proton-proton scattering. Blattnig et al. produced a parameterized 
spectral distribution by integrating E d 3 a/dp 3 (ref. 12). 

3. Results 

The Monte-Carlo programs JETSET and PYTHIA cannot be used in HZETRN due 
to excessive running time. Parametrizations of cross-section formulas are needed. The 
new idea is to make simple assumptions within the Lund Model framework to obtain 
analytic meson production cross-section formulas. Our purpose was not to produce any 
parameterization. To this end, all the fits in this work have parameters that are simply 
handpicked for the sake of illustrating the correctness of the qualitative aspects of cross 
sections predicted by string theory. Figures 1 and 2 clearly demonstrate the success of the 
fits. A major assumption in this paper is that a and b in equation (6) are not constants 
as proposed in the original Lund Model but are functions of rapidity. This assumption 
allows us to obtain the correct qualitative features of the cross sections without resorting 
to a Monte-Carlo simulation. At the present, the functions a and b are adjusted to fit the 
data. 

The exact forms of the parameters a and b (eqs. (9) and (10)) for various types of 
production are calculated from nonsupersymmetric string theory, such as the QCD string 
model (ref. 13). The concept of confinement by a string is quantified in terms of the 
minimal area law of the Wilson loop. QCD string models also include the gluonic degrees 
of freedom. The constants A and B in equation (13) and the constants c and k in equation 
(15) are also calculated from string theory. The (1 + 1) Lund Model cannot give angular 
dependence which is important for lowp z . Future work in string phenomenology hopefully 
can include the angular dependence of the cross sections near the threshold by extending 
the model to (3 + 1) or higher dimensions. Production formulas will also be extracted 
from JETSET and PYTHIA to be used in HZETRN. 
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Appendix A 
Kinematics 


The quark- antiquark pair of a meson is massless in the string model. In this case, the 
concept of the mass of a meson is associated with the mass of the string field and not the 
quarks. Massless quarks move at the speed of light. This result is not surprising because 
string theory predicts that the ends of open strings move at the speed of light either by 
imposing a Neumann boundary condition for open strings (ref. 14) or by solving the 
classical equations of motion (ref. 15). In the highly relativistic problems, the light-cone 
coordinates are the natural choice. In the (1 + 1) case, Lorentz transformations can be 
written as 


( E' \ = ( 7 -yu 

V P' ) V ~lfv 7 

( ? \ = ( 7 ~^ v 

y x' J y — qu 7 

where 7 is the Lorentz factor 


These equations imply 
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7 v/1 - v 2 


7 = cosh y p 
'yv = sinh y p 
v = tanh y p 


The energy and momentum of a particle can now be written as 



E = 7m = m cosh y p 
p = 'ymv = m sinh y p 

The boosted energy Ef, and the boosted one-momentum pb are given as 

E b = 7 (E-vp) 

= m (cosh y p cosh y — sinh y p sinh y ) 

= m cosh (y p - y) 

Pb = 7 (p~vE) 

= m (sinh y p cosh y — cosh y p sinh y) 

= m sinh(y p - y) 


(16) 

(17) 

(18) 


(19) 

( 20 ) 
( 21 ) 
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The relativistic velocity addition formula is simple in light-cone coordinates (ref. 5) with 


v' = tanh y' 

illustrating the additivity of rapidity, 


v ~ i’b 

1 — V b V 


tanh (y - y b ) 


y = y' + Vb 

This simplication motivates the definition of rapidity 
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Momenta along the light-cone can be defined as (ref. 5) 

p + = E + p = m cosh y p + m sinh y p = m e Vp 


P- = E — p = in cosh y p — m sinh y p = m e 


-yp 


With these definitions, boosts are simplified along the light cone (ref. 5), 

p' ± = m e ± ^ Vv ~ v ^ = p±e 1 y 

Light cone coordinates can also be defined in configuration space as 

in 

t = — cosh y„ 

K 

m 

x = — sinh y g 

K 


(26) 

(27) 

(28) 

(29) 

(30) 

(31) 

(32) 

(33) 


where k is the string constant and is used here to give t and x the correct dimension. 
A new subscript is adopted for y q in configuration space to distinguish y p in momentum 
space. These definitions are consistent with the requirement that v = dx/dt = tanhy 3 , as 
in equation (21). Similar results are obtained as in the momentum space case, 


m in m 
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m m m 
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Appendix B 
Lund Model 


The focus of this paper is primarily mesons. The force field between a quark-antiquark 
pair is presumed to be constant and confined to a flux-tube. This appendix is based on 
the work of reference 5, chapters 6 to 8. The equation of motion of any one member of 
the quark-antiquark pair acted upon by a constant force is (ref. 5) 


F = 


dp 

dt 


—k 


where k is the string constant. The solution of the equation is 


(37) 


P(t) = P() ~ Kt = Hi(t 0 - t ) 
From E 2 = p 2 + to 2 , E dE = pdp or 


(38) 


p dE 
E dp 

(39) 

is obtained. There is also 

dx dx 

p = 7 mv = 7 m — = E — 
dt dt 

(40) 

Equations (39) and (40) together yield 


dx p dE 

dt E dp 

(41) 

such that 

dE dE dp dt dp 

dx dp dt dx dt 

(42) 


Its solution gives the expected QGD flux tube result 

E = k(x o — x) (43) 

Combining equations (38) and (43) and m = 0, the relativistic equation of state is 

m 2 = E 2 — p 2 = k.[(x 0 — x ) 2 — (i 0 — t ) 2 ] = 0 (44) 


Let xq = to = 0, the equation of motion |®| = t is obtained. The maximum kinetic energy 
at x = 0 puts an upper bound on the displacement of the particle such that KX max = W . 
The path of the particle is illustrated in figure 3. The zig-zag motion of the particle is 
the reason for the name “ yoyo state.” The quark-antiquark pair of a meson are assumed 
to be massless. The mass square of the meson is proportional to the area of the rectangle 
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defined by the trajectories of the quark-antiquark pair. The period r of the yoyo motion 
is 

2Eo 
T = 

K 

where Eo is the maximum kinetic energy of each quark. In a boosted Lorentz frame, the 
period is transformed as 

t = t coshy (46) 

The breakup of a string occurs along a curve of constant proper time such that the 
process is Lorentz invariant. The Lund model assumes that the produced mesons are 
ranked, meaning that the production of the ?ith rank meson depends on the existence of 
rank n— 1, n — 2, ... ,1 mesons. A vertex V is a breakup point and the location k(x + , X-) 
where a quark-antiquark pair is produced. The breakup of a string is represented in 
figure 4. The produced particle closer to the edges are the faster moving ones correspond- 
ing to larger rapidities. 

Let p± o and p±j be momentum of the parent quark and the j th of the n-produced 
quarks moving along the x± light-cone coordinate, respectively. Then 



p± o = ^2p±j 

3 = 1 


and the momentum fraction at Vj is defined as 


z ±j 


P±j 
P± 0 


(47) 


(48) 


In order to simplify notations, Zj equals z +j unless specified otherwise, 
interval is 

ds 2 = dt 2 — dx 2 = dt 2 ( 1 — v 2 ) 

giving 

ds = dt\J 1 — v 2 = — = dr 
7 

where r is the proper time. The proper time is also x + x _ = (t — x)(t + x) 
and can be used as a dynamic variable such that 


The invariant 

(49) 

(50) 

= t 2 — X 2 = T 2 


1 1 2 
1 = K X+X- 


(51) 


Since T is proportional to the proper time square r 2 , its Lorentz invariant property makes 
it a good kinematic variable in the light-cone coordinates. Let W±i and W ± 2 be the 
kinetic energies along x± at vertices 1 and 2, respectively. The following identities can 
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be easily proven geometrically by calculating the areas of the rectangles F| = A, + As, 
r 2 = A 2 + A 3 , and m 2 as shown in figure 5: 


= (1 - Z -)W- 2 W +1 (52) 

r 2 = (1 - z+)W- 2 W +1 (53) 

m 2 = z-z + W-2W + i (54) 


where m is the mass of the produced meson. Equations (52) and (53) can be expressed 
with the help of equation (54) as 


„ 1 - z-) 

1 1 — 

z + z_ 

„ m 2 (l - z + ) 

1 2 — 

Z + Z_ 

Differentiating these equations gives 


dT j rn 2 

dz- z + z 2 

dT 2 rn 2 

dz + z\ z _ 


(55) 

(56) 


(57) 

(58) 


Let H(Ti) be a probability distribution such that H(I\) dT, dy, is the probability of 
a quark-antiquark pair being produced at the spacetime position Vi. From now on, the 
symbol V n also represents the breakup event that leads to the creation of the nth quark- 
antiquark pair along a surface of constant r. Let f(z + )dz + be the transition probability 
of obtaining V 2 given that V, occurs. The transition probability of V 2 via V\ is then 
H(Ti) dTi dy\ f(z + )dz + . Similarly the probability of V) via V 2 is H(T 2 ) dT 2 dy 2 /(z_)dz_. 
It is reasonable to assume that the probability of V\ via V 2 is equal to that of V 2 via V, 
such that 

iL(rx) dF, dy x f(z + )dz + = H( r 2 ) dT 2 dy 2 f{z„)dz_ (59) 

The Jacobian J in 

dT x dz + = J dr 2 dz _ (60) 


is 


J 


Mj. 

or. 

9r 2 

dz_ 

dz- 1 _ 

dz-\- 


1£Z 


£± 

z_ 


which can be easily computed by using equations (57) and (58) and <9Tx/<9T 2 
0. Equation (60) can now be reexpressed as 


dr 1 - 


dz. 


+ 


did 


dz_ 


z+ 


(61) 

dz + /dz_ = 

(62) 
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Since rapidity is additive ( i.e., y 2 = yi + A y) 


dy% = dy 2 (63) 

With equations (62) and (63), equation (59) is simplified as 

H(r 1 )z+f(z+) = H(r 2 )z-f(z-) (64) 

New definitions are now made to facilitate the solution of the equation 

h(T) = log if (T) (65) 

g(z) = log (z/(z)) (66) 

The new definitions transform equation (64) as 

M r i) + g(z+) = h(T 2 ) + g(z-) (67) 


Notice that 


d 2 g(z+) d 2 g(z _) 


0 


( 68 ) 


dz + dz- dz + dz- 

Differentiate equation (67) with respect to z^ and z_, eliminate terms with d 2 g/dz + dz_ 
using equation (68) and cancel out factors of <9T jdz on both sides of the equation to 
obtain 

5/i(ro , d 2 h(i\ ) _ dh(r 2 ) i r d 2 h( r 2 ) 


dT 1 


or 2 


or. 


or 2 


or equivalently 


d / dh 

dr v dr 


where b is a constant. The solution is 

h(T) = —bT + a In T + In C 

where C is a constant of integration. It yields the distribution 

H( T) = e h{r) = CT a e~ br 

Substituting equations (55) and (56) into equation (67) gives 
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where 512(2+) is 5(2+ ) of V\ changing to V2 and 521(2-) is g(z-) of V2 changing to V\. If 
a = a\ =a,2, the transition probability distribution is 

1 bmP 

•> ( 74 ) 

Z j 

where N is a constant of integration. Otherwise, the transition probability from V a to Vp 
is 

1 ( \ — Z \ a P — llr "' 2 

fap(zj) = N aP - -2- e ( 75 ) 

z i V z i ) 

where N a p is a constant of integration specific to the vertices V a and Vp. Let zoj be the 
jth rank momentum fraction scaled with respect to p + 0, then 


zo 1 = 21 ( 76 ) 

2b2 = 22(1—21) ( 77 ) 


The probability of two dependent events is the product of the probabilities of the two 
individual events. The existence of the rank - 2 hadron depends on that of the rank -1 
hadron according to the Lund Model so that joint probability of their mutual existence is 
the product of the two individual probabilities of I j and V2. With equations ( 75 ) to ( 77 ) 
the combined distribution of the rank-1 and -2 hadrons is 


f(zi)dz! f(z 2 )dZ2 


f(z m )dz 0l f 

N dz 0 i N dzo2 
201 202 

N dz 01 N dz 02 


202 


dzi 


02 


1 — 2oi / 1 — 2oi 

(l- 2 0 l) a (l- 


202 


201 


202 


. l—2oi 

(l-20l- 2 02 r e -^^) 


_ bm _ bm 

e z i z 2 


( 78 ) 


The geometrical identity Aj = bm 2 /zj is used in the last step. Generalizing the product 
of two vertices to that of n vertices, the differential probability for the production of n 
particles is easily seen as 


dP( 1 , . . . , n ) = (1 - z) a f[ e~ bA i 

J=i 


( 79 ) 


where 2 = YPj= 1 z o j ; let poj = zqj p + 0 and d 2 p = dp + dp _ , and use the identity 


J dCdBS(BC-D) 


dB 

~B 


to obtain 


dzoi dzo2 
2oi 2 q2 


cfpoi ^ 2 P02'5 + (Poi - m 2 ) ,5+ (Poi - m2 ) 


( 80 ) 

( 81 ) 
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With equation (81), equation (79) can be rewritten as (ref. 5) 

n 

dP(poi , . . . , Pon) = (1 - z) a N d 2 p 0j S + (plj - m 2 ) e~ bAj (82) 

i = i 

From simple geometry in figure 6, the kinetic energies of the quarks along the ± light 
cones are shown to be 


W- 

and the total kinetic energy square at V n is 


z Po+ 


E 


rri: 


i z OjPO+ 


s = W+W- = 


m 2 z 


3 = 1 ^ 

The total differential probability of n-particle production is (ref. 5) 
dP(z,s-, pou ■ ■ ■ ,POn) 


n o 

mz 


= dz5 z — Zqj ds S s — dP(poi , . . . , p$, 


3 = 1 


=1 Z 0j 


-S[l-y Uj J ds 5 - y '-Z- ) dP(poi, . . . ,P0n) 


■ i U 3 
J = 1 J 


where Uj = po+j/W n+ . Use the method 
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3 = 1 
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3 Prest ^ ^ POj 
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(83) 

(84) 
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(87) 
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where p rest = ( W n+ , W n J), to reorganize equation (86) as 


dP(z,s;poi,...,pon) 


ds - — (1 — z)~ e 


a ~~ br 5 


Prest 



x V d 2 poj d + (poj — m 2 ) e bArest 

3 = 1 


( 88 ) 


with the definition 

n 

Atotai = 'y ^ Aj = r + A rest (89) 

3 = 1 

The claim is that A to t a i is Lorentz invariant and is called the “Area Law.” (See ref. 7.) 
Equation (88) can be separated in the external and internal parts as in reference 5: 


dPext 


dPint 


day ( 1 - z) a e~ br , 

n 

13 N d^Poj 5 + (p 2 0j ~ in 2 ) S 
3=1 


Prest ^ ^ POj 


3 = 1 


g — bArest 


(90) 

(91) 


The external part contains kinematic variables s. z, and T. The internal part contains 
dynamic variables poj. Equations (90) and (91) are the final results. 
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X 


Figure 1. Experimental da/dx and string model results. Constants used in exponential function chosen to fit data 
and not parameterizations; p+p -*x° +X; P\ a \> = 360 GeV. 



Figure 2. Experimental E do/dpi and string model results. Constants used in exponential function chosen to fit data 
and not parameterizations; p + p — » Jt° + X; .P[ a b = 250 GeV. 
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Figure 3. Yoyo motion of quark-antiquark pair confined by linear potential. 



Figure 4. Breakup of quark-antiquark pair along surface of constant proper time x. Bold arrows represent velocities 
of produced mesons; breakup points labeled as vertices V\ to V n . 
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Figure 5 . Geometry of kinematics of two adjacent vertices V\ and V?, which are vertices or spacetime positions that 
also represent energy carried by string field. Quark moves along positive light cone and antiquark moves along 
negative light cone; antiquark of V\ combines with quark of Vi to produce meson of mass m\ W+\ is energy of V\ 
along x+; z + W + 1 is fraction of energy used to create quark from Vy, W _ 2 is energy of V2 along x_; z_W. 2 is fraction 
of energy used to create antiquark from Vp,A\,A2, A 3, and m- are areas of rectangles. From reference 5 . 


y's\ 



Figure 6. Geometry of kinematics of involving total area of spacetime diagram such that A tota i = T + A rest . Total 
energy square s = W+W . _ is represented by rectangle s; p± is energy of parent quark (antiquark) along ± light-cone 
coordinates; ZQ r p+ is fraction of energy used to create quark from V n ; m 2 /(zo n p + ) is energy used to create antiquark 
from V n -\. From reference 5 with minor modifications. 
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